Self-force of a point charge in the space-time of a symmetric wormhole 
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We consider the self-energy and the self-force for an electrically charged particle at rest in the 
wormhole space-time. We develop general approach and apply it to two specific profiles of the 
wormhole throat with singular and with smooth curvature. The self-force for these two profiles is 
found in manifest form; it is an attractive force. We also find an expression for the self-force in the 
case of arbitrary symmetric throat profile. Far from the throat the self-force is always attractive. 



I. INTRODUCTION 

Wormholes are topological handles linking different regions of the Universe or different universes. The activity 
in wormhole physics was first initiated by the classical paper by Einstein and Rosen [![, and later by Wheeler 
0. The latest growth of interest in wormholes was connected with the "time machine", introduced by Morris, 
Thorne and Yurtsever in Refs. 0, 0| • Their works led to a surge of activity in wormhole physics Q . The main 
and unsolved problem in wormhole physics is whether wormholes exist or not. The wormhole has to violate 
energy conditions and the source of the wormhole geometry should be exotic matter. One example [f| 0, HI of 
such exotic matter is quantum fluctuations, which may violate the energy conditions. Another possible sources 
are a scalar field with reversed sign of kinetic term (9| , and cosmic phantom energy [ToL [ll[ ■ The question of 
wormhole's stability is not simple and requires subtle calculations (see, for example [12|]). The problem arises 
because a wormhole needs some amount of exotic matter which violates energy conditions and has unusual 
properties. Recently it was observed that some observational features of black holes can be closely mimicked by 
spherically symmetric static wormholes [l3j having no event horizon. Some astronomical observations indicate 
possible existence of black holes (see, for example, Ref. [H|]). It is therefore important to consider possible 
astronomical evidences of the wormholes. There is an interesting observation of quasar Q0957+561 which 
shows the existence of a compact object without an event horizon. Some aspects of wormholes' astrophysics 
were considered in Ref. [l6| , where it was noted that this massive compact object may correspond to a wormhole 
of macroscopic size with strong magnetic field. A matter may go in and come back out of the wormhole's throat. 

In the framework of general relativity there exists a specific interaction of particles with gravitating objects - 
the gravitationally induced self-interaction force which may have considerable effect on the wormholes' physics. 
It is well-known [17| that in a curved background alongside with the standard Abraham-Lorenz-Dirac self-force 
there exists a specific force acting on a charged particle. This force is the manifestation of non-local essence of 
the electromagnetic field. It was considered in details in some specific space-times (see Refs. [HI, Gil for review). 
For example, in the case of the straight cosmic string space-time the self-force appears to be the only form of 
interaction between the particle and the string. Cosmic string has no Newtonian potential but nevertheless a 
massless charged particle is repelled by the string [2(j, whereas massive uncharged particle is attracted by the 
string due to the self- force [2l[ . The non-trivial internal structure of the string does not change this conclusion 
[22| . The potential barrier appears which prevents the charged particle from penetrating into the string. For 
GUT cosmic strings the potential barrier is ~ l0 5 Gev. The wormhole is an example of the space-time with 
non-trivial topology. The consideration of the Casimir effect for a sphere that surrounds the wormhole's throat 
demonstrates an unusual behavior of the Casimir force [23[ - it may change its sign depending on the radius 
of the sphere. It is expected that the self-force in the wormhole space-time will show an unusual behaviour 
too. As we will show in next sections the self-force indeed has an unusual behaviour - the charged particle is 
attracted by the wormhole instead of expected repulsive action. This result is valid for arbitrary profile of the 
throat. Therefore charged particles are likely to gather at the wormhole's throat. 

The organization of this paper is as follows. In Sec. [IT] we briefly discuss the background under consideration 
and describe all the geometrical characteristics we need. In Sec. IHII we develop a general approach to calculation 
of the self- force in this background for arbitrary profile of throat and apply it to two specific kinds of profiles. 
We obtain simple formulae for the self-force in these cases. However, this approach is valid only for simple 
enough throat profiles, and therefore in Sec. IIVI we develop an alternative method which is appropriate for an 
arbitrary symmetric throat profile. We obtain general formulae and find general expression for the self-force far 
from the wormhole's throat. We discuss our results in Sec. |Vj Throughout this paper we use units c = G = 1. 
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II. THE BACKGROUND 

Let us consider an asymptotically flat wormhole space-time. We choose the line element of this space-time in 
the following form 

ds 2 = -dt 2 + dp 2 +r 2 {p)(d0 2 + sin 2 9dip 2 ), (1) 

where t, p £ R and 9 £ [0, ir], ip £ [0, 2ir}. Profile of the wormhole throat is described by the function r(p). This 
function must have a minimum at p = 0, and the minimal value at p — corresponds to the radius, a, of the 
wormhole throat, 

r(0) = a, f(0) = 0, 

where an over dot denotes the derivative with respect to the radial coordinate p. Space-time is naturally divided 
into two parts in accordance with the sign of p. We shall label the part of the space-time with positive (negative) 
p and the functions on this part with the sign "+" ("—"). 

The space-time possesses non-zero curvature. The scalar curvature is given by 

2(2rf + f 2 - 1) 

it — 2 • 

Far from the wormhole throat the space-time becomes Minkowskian, 

Kp)Ip-±°o =±p. (2) 

Various kinds of throat profiles have been already considered in another context [7| . The simplest model of 
a wormhole is that with an infinitely short throat [7[ , 

r = a+\p\. (3) 

The space-time is flat everywhere except for the throat, p = 0, where the curvature has delta-like form, 

r = -s 6 M. 

a 

Another wormhole space-time that is characterized by the throat profile 



r = Va 2 + P 2 (4) 

is free of curvature singularities: 

R 2a 2 

(a 2 + p 2 ) 2 

The wormholes with the following profiles of throat 

r — p coth — + a — b, 
b 

i P 

r = ptanh — + a. 
b 

have a throat whose length may be described using a parameter b. The point is that for p > b the space-time 
becomes Minkowskian exponentially fast. 

III. SELF-ACTION IN THE WORMHOLE SPACE-TIME 

Let us consider a charged particle at rest in the point p',0', <p' in the space-time with metric (Q}. The Maxwell 
equation for zero component of the potential reads 

AAQ = W(p-p') 6(0 - 0')6(<p - <p') 
r 2 (p) sin 9 

where A = g kl S7k^7i- Due to static character of the background we set other components of the vector potential 
to be zero. It is obvious that A° = 47reG(x;x'), where the three-dimensional Green's function G obeys the 
following equation 

Art >\- Kp-p')W-0')K<P-<f f ) 

AG(X ' X) ~~^(^ ri^e ' 
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Due to spherical symmetry we may extract the angular dependence (denoting succinctly £1 = (9, (p)) 



oc I 



G(x;x')=]T Yim{WLW)m{p, P ') 



1=0 m=-l 

and introduce the radial Green's function gi subject to the equation 



2r' . 1(1 + 1) 



9i H 9i 

r 



-91 



s(p-p') 

r 2 (p) 



We represent the solution of this equation in the following form 

gi = 6{p - p')^x(p)Mp') + 8(P' - p)*i(p')*2(p), 
where functions ^ are the solutions of the corresponding homogeneous equation 

1(1 + 1) 



T 2r ' T 

* H * 

r 



= 0, 



satisfying the boundary conditions 



lim \l>i = 0, lim * 2 

p — >+00 p — >+oc 



If one substitutes ([6]) to (J5]) the condition for the Wronskian emerges: 



(5) 



(6) 



(7) 



(8) 



(9) 



We consider the radial equation in domains p > and p < and obtain a pair of independent solutions 
(j) 1 , <fi 2 for each of the domains separately. We do not need to consider two domains if it is possible to construct 
solutions that are C 1 -smooth over all space. However, this is not the case for many situations. For the two 
kinds of the throat profile considered below we may easily construct solutions for the two domains separately 
(but not for all space). After that a procedure developed here allows to construct C 1 -smooth solution over all 
space. Due to condition ([2]) we obtain asymptotically 



j 



<P Ip^+oo = P 

The Wronskian of these solutions has the following form 



-i-i 



w(4> l ± Al) 



A± 
r 2 (p) 



(10) 



with some constants A±. Let us consider two different solutions: 

*i = 

*2 = 



aUl+pU 2 ,, p>0 



,J + 

alcfrl + f3\<t> 2 _, p < ' 

a\<l>l + 0l<fi., p>0 
a 2 _<$>_ + P_(t> 2 _, p < ' 

The Wronskian condition (fT0|) implies the constraint on the coefficients: 



The general solution of the matching conditions 



a i = a_ 



W(<f>_, 


2 +) 


W(4\, 


</> 2 +) 


w(4>\, 




w(4>\A\) 



*+(0) = *_(0),*+(0) = #_(()), has the following form 



W(<t>\,4? + ) 

W(4l,^) 



w(<j>\^+) 
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or vice versa 



0- = —a+ 



w{4>\,4> 2 -) , a w{4>\^ 2 _) 
+ p+ 



w{4> l _,(t? 



■fa 



w{<fii,<jrL) 



To satisfy the conditions ([5]) we consider two specific solutions that emerge when we set a\ = and o?_ = 0. 
Corresponding solutions have the following form: 

vT, - J 44' P > ° 

1 S aL^L +/3i0 2 _ ! p < ' 



ai^i+tf^, p>0 



where 



and 



„i fl i W^ 2 -) 



2 

a + — P- 



,2 _. , „2 WtflM 



w(<^>4) 



The Wronskian condition reads 



Taking into account the above relations we obtain the radial Green's function in the following form 
1. p > p' > 



ffPW) = --^#V)#.(p) + 



i 



i WW.) 



.4 



2. < p < p' 



A + W{4>\A 2 -) 
g\ 2 \p,p')=g\ x \p',p) 



3. p < p' and p' > 0, p < 



i w&x^D 



4. p > p' and p' < 0, p > 



5. p' < p < 



~ A + W(4> 2 + ,<f> 2 _ 
gl 4 \p,p')=gl 3 \p\p') 
'\ 1 a?( 'W/ \_l 1 ^-»4) 



<£(pV-(p) 



6. p < p' < 



ffl (8 W)=fffV,P) 
where the constants A± may be found from the relation (fT0|) : 

at arbitrary point p. 



(11a) 
(lib) 
(11c) 

(lid) 
(lie) 

(llf) 
(12) 
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Let us consider in detail the simple case of the symmetric throat profile: r{—p) = r(p). In this case we may 



choose (j>]l 2 (p) = ^Kl 2 (—p) and hence (0) = </>i' 2 (0) and ^^ 2 (0)' 
account these formulas we obtain 
1. p > p' > 



^ 2 (0)' and A + = Taking into 



1 

X 



i w + 



2. < p < p' 



(2) 



(p,p') =g i i 1 \p',p) 



4>\{p')4>1{p) 



3. p < p' and p' > 0, p < 



4. p > p' and p' < 0, p > 



sfW) = - 



1 I4^,0 2 +) 



2 + (p')^»+(-p) 



5. pf < p < 



6. p < p' < 



(4) 



iP,P) 



(3) 



(/>V) 



■p,-pO 



2/12/2 + 2/12/2- Thus we have to write out in manifest form only and <7, . 



(5)/ 



(13a) 



(13b) 



(13c) 



(13d) 



(13e) 



(13f) 



Here W+(y 1 ,y 2 ) 

The Green's function will not give a finite expression for the self-force. The origin of this divergence is 
the electromagnetic self-energy of point particle. There exist some approaches to obtain finite result. First 
simple way is to consider total mass as observed finite mass plus infinite electromagnetic contribution. Usually 
this subtraction is called classical " renormalization" because there is no Planck constant in the divergent term. 
Dirac [2(| suggested to consider radiative Green's function to calculate the self-force. Since the radiative Green's 
function is the difference between retarded and advanced Green's functions, singular contribution cancels and 
we obtain a finite result. There is also axiomatic approach suggested by Quinn and Wald [27| . In the framework 
of this approach we obtain finite expression by using a "comparison" axiom. This approach was used in Refs. 
[28l . |29| for a specific space-time. 

We will use general approach to renormalization in curved space-time [2~i| which means subtraction of the 
first terms from DeWitt-Schwinger asymptotic expansion of a Green's function. In general there are two kinds 
of divergences in this expansion, namely, pole and logarithmic ones plj |. In three-dimensional case which we 
are interested in there is only pole divergence, while the logarithmic term is absent. The singular part of the 
Green's function, which must be subtracted, has the following form (in 3D case) 



G 8 



1 A 1 / 2 
4tt V2^' 



where a is half of the square of geodesic distance and A is DeWitt-Morrett determinant. If we take coincidence 
limit for angular variables then these quantities are easily calculated using the metric (fTJ) : a = (p — p') 2 /2 and 
A = 1. Therefore to carry out renormalization we have to subtract from the Green's function its singular part, 
which has the following form: 



1 



1 



4ir \p - 



This is, in fact, the Green's function in Minkowski space-time. This approach was used many times in different 
curved backgrounds (e.g. see [13]). Now we are in a position to consider first some specific cases in detail and 
after that to proceed to the general profile of the throat. 



G 



A. Profile r — a + \p\ 



Two linearly independent solutions, (j) 1 , </> 2 , are given by 

r(±p) 1 = (a±p) 1 , 

a 2l+1 r{±p)- l - 1 = a 2l+1 {a±p)- l -\ 
with the Wronskian 



Therefore the solutions we need have the following form 



r 2 (p) ' 











= 


2(*+l) <1 






I 2Z + 1 V- 








*2 




2/+1 <?+ 



J_ ^2 



p > 
p<0 ' 

p > 
p < 



From Eqs. @ and (fTS")) we obtain an additional constraint 

1 



a 2l + 1 2(l + 1) 



Thus we have found the radial part of the Green's function in manifest form 



(2l + l)g ( l 1 \p,p') 
(2Z + l).gPW) 



,22 + 1 



rV) 1 

r l + x {p) 2(Z + l)r i + 1 ( j o)r i + 1 ( j o')' 



2Z + 1 



,2/+! 



2(/ + l)r i + 1 (p)r i + 1 ( ( o')' 
To calculate the full Green's function and the potential we turn to the relation 

l 

47T 



£ Y lm (Q)Y l l l (n') = ^P l (cos 1 ), 



m — — l 



by using which we obtain 



oo 



1=0 



(14a) 
(14b) 



(15) 



(16a) 
(16b) 



Here cos 7 = cos 8 cos 8' + sin 8 sin 0' cos(<^ — ip') is a cosine of an angle between two points on a sphere. Now we 
use two series expansions: 



1 



fc=0 

t k+l 



y/l - 2tx + t 2 



hi 



fe=0 



1 + 



2t 



1 - t + Vt 2 - 2tx + 1 



which allow us to obtain the Green's function in manifest form 



where t 



4irG {1) {x;x') = 

4irG {3) (x;x') = 
f 

r(p)r(p') ■ 



1 



y / r(p) 2 - 2r{p)r(p') cos 7 + r{p') 2 2a 
t 1 



2t 



l-t+ ^t 2 - 2tcos7 + 1 



iVi 2 - 2icos7 + 1 2a 



hi 



1 



2t 



1 -i + v/t 2 -2tcos7 + l 



FIG. 1: Several surfaces of constant potential are shown. The charge e = 1 is at the point x — 2,y = z — 0. The small 
sphere is the throat of the wormhole, p — 0. We observe that some surfaces go under the throat to another universe. 



Let us consider the case when the observation point, p, is far from both the wormhole throat and from the 
position of a particle, p' , that is we set p — > oo. We consider p' > 0. In this limit we have 

p > : 4ttG(x; x') = - - ° + 0{p- 2 ), (17a) 
p 2p{a + p') 

p < : 4ttG(x; x') = - ° + 0(p- 2 ). (17b) 
2p(a + p') 

From the first of these expressions we observe that the potential of the particle contains an additional term with 
the same p dependence alongside with the standard Coulomb part. The same problem appears for a particle in 
the Schwarzschild space-time [3l[ . To solve this problem it was suggested to add a solution of the corresponding 
homogeneous equation to the solution with a "bad" asymptotical behaviour (32J. However, in our case there 
is no need for such modification: the potential (jTTJ) correctly describes real potential of the particle at rest. 
Indeed, let us consider the Gauss theorem in the space-time of the wormhole. We place the particle at the point 
p' and enclose this charge by two spheres, namely, S± : p = R and 5*2 : p = —R. The Gauss theorem reads 



4vre = JJ EndS* + JJ EndS 1 , 



Si s 2 

where E = — VAo and n = (±1,0,0) for S\ and S*2 correspondingly. Therefore we have 

47re = - ffd p A Q \ p=R dS + JJ d p A \ p= _ R dS. (18) 



Si S2 



Taking into account formulae (|17|) it is easy to show that the Gauss theorem is satisfied and we do not need to 
add a solution of the homogeneous equation to correct the potential. Thus the potential has "bad" behaviour 
at wide separations of the observation point and the particle 



A Q = AneG{p; p 1 ) w - 
P 



1 



2{p> + a) 



instead of the expected expression e/p. In fact, we observe the particle with charge e(l — 2 ( P ^+a) ) ^ ar f rom 
our genuine charge e. The explanation is very easy - some electric field lines have gone through the throat to 
another, invisible domain of space-time and for correct formulation of the Gauss theorem we must take these 
field lines into account, too. Hence the second sphere S2 emerges. The Gauss theorem relates the charge of the 
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particle to the electric field flux density. In the Fig. [T]the surfaces of the constant potential Ao of the particle 
are plotted. 

Now let us restrict ourselves to the + part of space-time and consider the situation when f) = Of. We have 
the following expression for the Green's function 

III " 2 
G{p-p') = -- - + _ln[l- 



4-7T |p — p'\ 87m r(p)r(p') 
Therefore, after renormalization we obtain 

G ren (p; P ') = ^-Ml 



87ra r(p)r(p') 
The potential at the position of the particle reads 



* = £ln[l- 



2 J 



and the self-potential is 



Some limiting cases are: 



2a (a + p) 



U=— Hi--, -]. (19) 

4a L {a + p) 2 

^U=o = -S + ^» (2°) 



U\ a ^0 = 



Ap 2 2p 3 ' 

2 

ae 
4(a + p) s 



The self-force 



has the radial component only 



F = -Vt/ 



F" = -d p U 



ae 2 



2r(p) 3 1 



Kp) 2 



The self-force is always attractive, it turns into infinity at the throat and goes down monotonically to zero 
as p — > oo. We may compare this expression with its analog for Schwarzschild space-time with Schwarzschild 
radius r= = a: 



l- a - 
2 r 3 y r 2 ' 

Important observations are: 

1) The self-force in the wormhole space-time has an opposite sign - it is attractive. 

2) Far from the wormhole throat and from the black hole we have the same results but with opposite signs 

F" = -— 

wh 2f)3 . 

3) At the Schwarzschild radius r s — a the self-force equals zero, whereas at the wormhole throat it tends 
to infinity. The latter discrepancy originates in the selected throat profile function that leads to the curvature 
singularity at the throat. 
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B. Profile r = a 2 + p 2 

The radial equation ((TJ) has the following two linearly independent solutions 

^ = cfPiiz), (21a) 

4>\ = c+Qi(z), (21b) 

with Wronskian 

W{<j> l + ^l)=i44^. (22) 

Here Pi and Qi are the Legendre polynomials of the first and second kind, and z — ip/a. In the particular case 
when 4 = (2i-i)\\ > c 2 ~ TZ^yrrrfj the solutions are real functions and far enough from the wormhole we have 

,2 I 21 + 1 -l-l 

= a p 

as in the above example. We should note that solutions given by Eq. (j2Tj) are not C 1 -smooth over all space. 
The point is that the above functions have the following symmetry properties [331 ] 

Pii-z) = (-l) l Pi(z), Qi(-z) = (-l)' +1 Q ; (z). 

These conditions imply that the functions are not C^-smooth at the throat p = 0. However, the procedure 

developed above yields C 1 -smooth solution provided that we have solutions in domains p > and p < 0. For 
p < we choose 

<t>- = c-pi-z), (23a) 

<t>- = cZQt-z), (23b) 

with Wronskian 

W(^,<j> 2 _) =-ic[<5± (24) 

To simplify the expressions we take c± = cf , = '. Final result does not depend on the specific choice 
because the formulae involve only the ratio of function values at points p and p = 0. 

To proceed further we use the formulae (we take into account the limit as p — > along the imaginary axis 
from above) 

P(0) = ; n ,P{(0)= ^ 



r(|-i)r(i + i)' r(| + i)r(-i)' 

2 e r(i + i)' WUj Vne ra + iy 



With the help of these relations it is easy to show that 

W+{cj>\,<j>\ 



I c{ 

7TC+ 

- r — - - , _i ,/ - ill' 



(25) 



= (-l) i+1 -^, (26) 



^ C 2 



A± = iicfcja. (27) 
Now we use the formulae (fT5j) and (f2"4"| to obtain the radial Green's function 

5i (1 W') - -Pi{z')Qi{z) + —Qi{z)Qi{z'), (28) 
a ira 

gf\p,p') = —{-l) l+1 Qi{-z)Qi{z'). 



ira 

Therefore 



G {1) (P,P') = j- f> + 1) \iPi(z')Qi(z) + -Qii^Q^z)] 
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and 



1=0 



1 00 

where z = ip/a and z' — ip'/a. Now we use the Heine formula [33] 

jyil + l)P l {z')Q l {z) = — 

and obtain 

4ttGW( P ,p') = — jT(2Z + l)(-l) i+1 Q,(z')Qz(-*) 

IT/1 fc * 



1 1 °° 

i , + -32i + i)QiMOiW, 

|p — p i 7ra z — ' 



To find these series in the closed form we use the integral representation for the Legendre function of the 
second kind 13311 



Qi(z) 



i f 1 W) 



z - t 



d,t. 



which yields 



J2{2l + l)Qi(ix')Q l (ix) = - 

1=0 

oo 

^(2Z + l)(-l)' +1 Q^0Qz(-^) = - 



1=0 



arctan X — arctan x' 
x — x 

— arctan x + arctan x' — 7r 
—x + x' 



Therefore 



AnG^(p,p') 



1 



1 arctan - — arctan B- 



\p-p'\ 7T p-p' 
1 — arctan £ + arctan £- — 7r 

7T ~P + P' 



Far from the wormholc throat we have 



lim 4ttG (1 < 2) (p,p') = - + - 

p^+oo p p 



1 1 p' 

1 — arctan — 

2 7r a 



lim 47rG (3) (p,p') = - 



1 1 p' 

1 — arctan — 

2 7r a 



hence the potential far from the charge kept at the point p' has the following form 

, n e e 

A °(i) = - + - 
p p 

e 



A 



1 1 p' 

1 — arctan — 

2 7r a 



1 1 p' 
— I — arctan — 

2 7r a 



and again it obeys the Gauss theorem (|18p. The renormalized Green's function is given by 

1 °° 

47rG(p,p) ren = _^(2/ + l)Q^) 2 



Z = 

1 a 



7r p 2 + a 2 
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Therefore we obtain the self-potential 



and the self-force 



U = 



2ir p 2 + a 2 



9 ir (p 2 + a 2 ) 2 



(29) 



(30) 



As expected the self- force is everywhere finite and equals zero at the throat. Far from the wormhole we have 



pp 



e 2 a 



7T p° 



Thus the self-force is always attractive. It has maximum value at distance p* = a/\/3 with magnitude F^ aa , = 
3v / 3e 2 /167ra 2 . The plots of the potential and the self-force are shown in the Fig. [2l 



U a F p , 




FIG. 2: Thin line is plot of potential (|29|l and thick line is a plot of the self- force (I30|) . The force has an extremum at 
point p* — a/ -\/3. 



IV. GENERAL CASE 



In this section we outline our approach to generalization of the results of the previous sections. From now 
on the profile function r(p) will be considered to be arbitrary symmetric function (the condition of asymptotic 
flatness given by Eq. (J2]) is understood). 

Let us perform the WKB analysis of the defining equation for radial Green's functions 

J + ^_J2 + 1>, = . (31) 



(32) 



To accomplish this we represent the solution of this equation in the form 

and arrive at the following equation for 5": 

§ + S 2 + yS - ^ ~ 2 1/4 = 0, (33) 
where v = I + 1/2. Next step is to expand S in the following power series in v [30| : 

oo 

S = v ~ ns n- (34) 
n=-l 

One comment is in order. Equation ()33|1 is valid for arbitrary I > 0, whereas the expansion over \ jv ((34j) 
is fulfilled for I > only. For I — we have 1/v = 2 and the expansion is divergent. Nevertheless, we may 
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easily find gi for / = without turning to the series expansion. Indeed, for I = the equation (|3Tj) simplifies 
considerably (ip stands here for the zero mode only): 

2f 

(p -\ ip = 0. 

r 

This equation has two independent solutions: 

<p\ = C\, p\ = C 2 f \dp + C 3 . 

J On ' 



We take C% = 1 in order to be in agreement with (|14p but this is insufficient. This solution corresponds to the 
solution ~ p l for I — 0. The second solution must tend to zero far from the throat and it corresponds to the 
solution ~ p 11 for I = 0. Therefore we may write 



poo 

lim <p 2 + (p) = C 2 / —dp + C 3 = 



hence 



-Co 



7 dp. 



Thus we obtain the second solution in the form 

IP%{P) = -C2 



-dp. 



We set C?2 = —a to be in agreement with (TT4")) . To summarize, we have 



1, ¥>+ 



7 dp, 



(35) 



with Wronskian W(</?j_, 
function 



-a/r 2 . These expressions allow us to obtain zero component of the Green's 

y + (p') 



i 

2a'' 



Now that we have separated the I — case, it is possible to proceed with I > contributions, using the series 
expansion. Substituting ([34]) to |33|) yields the following set of recurring expressions for the functions S n : 



S-i 
Sq 

Si 

Sm+1 



1 

±-, 

r 



S-i r 
1 



25_i 
1 

25-1 



2f 



1 



S + S + —So + — 



2r 



m = 1,2,. 



The sign "+" in the first equation corresponds to the solution 1 and sign "— " to the solution <fi 2 which tends 
to zero as p — > oo. Following are the manifest expressions for derivatives of the first several functions S n : 



S-i 


1 

= ±- 

r 




So 


r 




Si 


1 

= H 


r 2 1 
r + 27 ~ 27 


S 2 


i 

~8 


rr (3) + 2rr 



l f ■■ 1 -2v 

= -g(rr + -r ) , 
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16 
1 

32 



r 2 r ( 4 ) -|- 4 r) i r ( 3 ) _|_ ^ r f 2 _|_ + _f: _ !L + T _ 

2 2 2 8r 4r 8r 



r 3 r (5) + 7 r 2^ r (4) + 6r 2 fr (3) + 9rfV 3) + rr (3) + 4rrr 2 + 2ff 

= — (4rV 4) + 16r 2 rr (3) + 4r 2 f 2 + 4rr 2 f + 4rf - r 4 + 2r 2 ) 
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Note that all £271 come with only one sign, whereas the signs of S2n+i vary according to the choice between <f> 
and (j> 2 . Another important point is that all S^n appear to be full derivatives, in contrast to all S^n+i- These 
two points will be essential in the subsequent analysis. 

A simple verification of the approach can be done at this point. Taking r = \J p 2 + a 2 we can evaluate 
the expansion and then compare (|3"2"1) with the exact solutions (|2"Tj) . Numerical calculation reveals fast 
convergence of the series to the exact solutions. The series that results when all S^n+i are taken with the "— " 
sign coincides with the solution i n+1 Qi(z) (|21b[) . while the series with the "+" sign corresponds to the following 
linear combination of (|21a[) and (|21b[) : (— i) n (Pi(z) — -Qi(z)). We would like to note that the result of the 
previous section would be the same if we had used this linear combination instead of (—i) n Pi(z). 

Let us introduce a pair of indices to help us distinguish between different solutions. From now on S 1 will 
denote the whole series (jMJ), in which all S^n+i are taken with the "+" sign, while S 2 will correspond to the "— " 
sign. These are essentially the two different solutions that we previously called tfi 1 and <fi 2 . We will additionally 
mark the solutions with "±" according to the sign of p. The Wronskian condition (TIT))) is therefore cast into the 
form 

,S±\ P )+S ±2 (p) A ± 1 



r 2 (p) S ±2 {p) - S±\p)' 



and evaluating it at p = we get 



a 2 5± 2 (0)-5 ±1 (0)" 
Functions S^n with even indices have the following structure 

n 

s 2 n(p) = X>* +1 r( 2 * +1) ; 

each term necessarily contains odd order derivatives of r with respect to p. Hence if the profile function r(p) is 
symmetric, r{—p) = r(p), then S^^O) = 0. 

By using the above formulae we may represent the Green's function (|13ap in the following form 



i S+ 2 (p)-S+ 2 (0)+S+V)-S +1 (0) 

gr(p,p') = — j = = (36) 

1 e S+ 2 (p)-S+ 2 (0) + S+V)-S+ 2 (0) 5-2(0) _ g+l^ 

+ a 2 S+ 2 (0) - S+i(0) S- 2 (0)-S+ 2 (0)' 

Due to symmetry of the profile function, r(—p) = r(p), we may replace S~ 2 (0) = — 5 +2 (0). From the above 
chain we obtain in manifest form: 



Si'j^-SiftO) = ± 



p dp 
r{ P y 



^ 2 (p)~^o' 2 (0) = iln^. 
Therefore the radial Green's function reads (I > 0) 

1 e ~ u K'^ e 12™ =1 v- n {si 2 {p)+st 1 (p')) 



9i(P,p') 



^vkpW) e: = o^"«-i(o)) 

^ V^pHp 7 ) Er=o^ 2 "(«4 + n-i(o)) E,T=o^" +1 (-«^ 2 i(o))' 



14 



The second part here gives zero contribution because S 1 ^ 1 contains only odd order derivatives of r, which are 
zero at p = 0. This fact may be verified straightforwardly by performing calculations for the \/p 2 + a 2 profile 
function. In this case 

S +2 (0)+5 +1 (0) =0. 

We must take into account that S +2 corresponds to Pi(z) and S +1 corresponds to Pi(z) — -Qi(z). This issue 
was already addressed above. 

Let us now perform summation over I in order to obtain full Green's function. We have 



1=0 



1 00 1 



where 



1=1 1=0 



dp 



P ' r(p)' 

To obtain the above formula we have permuted the summations over k and over The functions //- are expressed 
in terms of the function <E> [33| 



/ fc (6) = e ->$(e- 6 ,A l -; 



First two functions in manifest form are 



1 1 



2 smh(|) 

Mb) = lni±^4-2e-^. 
1 — e 2 

To calculate the self- force we need the coincidence limit p' —> p. In this case only first two function are 
divergent 

fo(b) 1 1 n , is 

--- + 0(p-p), 



VhpMT) p-pf 

fi(b) 1, p-pf 2 

In— \-0(p- p ), 



r 



~ -(H{k,~) + U(p- p), 



y/ripW) r 2' 

where is the Hurwitz zeta- function [331 ] - 

The functions jk can not be found in closed form for arbitrary index fc, however, each function may be found 
from the expansion over v. In manifest form we obtain the following expressions for these functions: 

Jo(p,p') = 1, 

•/ /\ f p -l + r 2 + 2rr -1 + r 2 + 2rr f .. 2 

3i{p,p) = - / s dp= (p - p') + 0((p - p'y), 



jpi 8r r 8r 

. , -l + r 2 + 2rr 3 2 
3i[p,p) = g = gflir , 

j 4 (p,p) = — [3 + 3f 4 -12rf-4r 2 f 2 -2f 2 (3 + 2rf)-32r 2 fr (3) - 8rV 4) ], 
128 

J 6 (p,p) = Y^-[5-5f 6 -30rf-20r 2 f 2 +3f 4 (5 + 6rf)-192r 2 r 3 r (3) 
- 8r 3 (f 3 + 5r (4) ) - f 2 (15 + 20rr + \2r 2 r 2 + 456rV 4) ) 
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8r 4( llr (3)2 + 20f r (4) ) - 16r 2 r(10(l + 3rr)r {3) + llrV 5) ) - 16r 5 r 



(6)1 



In the above expressions a\ stands for the first heat kernel coefficient (see review |34j|). 
All terms of the form /2fe+ij2fc+i vanish in the coincidence limit. Therefore we obtain 

00 2u 1 1 11 °° 3 11 

f— ' 47T 47r — p' r r f—* 2 a 2a 



2=0 



fc=i 



V\{P) 

V+(o) 



After regularization we arrive at the formula 



u(p) = y 



fe=i 



This expression is exact and we may use it for arbitrary throat profile. 



V 2 + (o) 



(37) 



-0.06 



-0 .08 




FIG. 3: Thin line is the plot of the exact formula for the self- force given by Eq. (|30[) . Middle thickness line is the plot 
of the self-force calculated by the above formula up to the first term k — 1. Thick line is the same formula up to k = 2. 



Fig. [3] shows the exact form of the self-force for the second kind of profile of the throat and the result obtained 
from the above formula. We observe that even the first term of the series gives good approximation. In fact 
both expressions are very close for small and large distances from the wormhole throat, there is small deviation 
close to the extremum only. For this reason we may write out an approximate expression for self-potential in 
which there is contribution from first heat kernel coefficient 



U(P) « - 



1 , /„ 3 3 2 1 a 1 a y+(p) 



(38) 



where Ch(2, f) = ^--4si 0.9348. 

We may obtain general conclusion from (|37|) . Let us consider large distances from the throat and suppose 
that the profile function has the following asymptotic series expansion: 



P- 



VnP 



Each coefficient j2k (p, p) has the same asymptotic form 

32k(fi,p) = -{^ + 0(p~ 3 )- 

The relation 



00 ^ 4 



fe=i 



2' 3 



yields 



1 q 

- J2 (H(2k, -)j 2k (p, p) = + 0(p- 4 ) 



fe=i 



1G 



Putting all expressions together we obtain the following asymptotic expansion for the renormalized Green's 
function 



and for the self-potential 



G ren ( P ;p) 



U = — 



1 a 



1 avo 



8irp 2 ip 2 + (0) 4np 3 ip 2 + (0) 



0(p~ 



e av 



+ o( P - 4 ). 



V^( ) 2/^(0) 

Let us compare this result with the cases considered above. For the profile r = a + \p\ we have v° = a and 

adp 



o (a + Pf 



= 1. 



Therefore 



U 



9 9 9 



Ap 2 2p 3 

This expression coincides with the result that follows from the exact formula (f2T))) . For the smooth throat profile 
y/ p 2 + a 2 we have v° = and 



^+(0) = 



adp 7r 



o a 2 + p 2 2 ■ 



Therefore 



U 



2np' 



+ 0(p~ 4 ) 



in agreement with Eq. (|29[) . 

Thus we may conclude that far from the wormhole throat we have the following expression for the self-potential 



U = 



4p 2 <p\ (0) Ap 2 



dp 
r 2 (p) 



Note that it is always negative, hence the self-force is an attractive force. All information about the specific 
throat profile is encoded in the factor 

dp 
r 2 (p)' 



V. CONCLUSION 



The aim of the paper was to calculate the self-force acting on a point charge at rest in the space-time of a 
wormhole. This calculation has relevance to many problems. First of all the space-time of a wormhole is an 
example of topologically non-trivial space-time and it is interesting to clear up the role of non-trivial topology 
of the wormhole space-time. As far as we know this is the first calculation of this kind. The second point is 
that the self-force may play a serious role in wormhole physics. 

We have developed standard approach to calculation of the self-force. Following it we have obtained the 
self-force explicitly for two profiles of the wormholc's throat. The first shape of the throat is characterized by 
a curvature singularity - the space-time is everywhere flat except for the throat. The self-force is everywhere 
attractive and singular at the throat. The singularity originates in the accepted model of the throat. Far 
from the wormhole's throat the self-force has the same magnitude as in the Schwarzschild space-time with the 
Schwarzschild radius equal to the throat radius, but the sign of the force is opposite. In the second example we 
have carried out the calculations for a smooth throat profile. The self-force is zero at the throat, as expected 
due to symmetry of the space-time with respect to the throat. In both examples we observe that the self-force 
is attractive, irrespective of the particle position. We suppose that this is inherent to the wormhole with an 
arbitrary throat profile and is the manifestation of non-trivial topological structure of space-time. The manifest 
calculations of self-force for arbitrary throat profile developed above show that far from the throat the self-force 
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is always attractive and falls down as third power of distance. The magnitude of this force depends on the 
throat profile. 

Let us now speculate about the result obtained. Charged particle in the wormhole space-time is always 
attracted towards the wormhole's throat. Starting from rest at infinity it will reach the throat with kinetic 
energy (for the second kind of throat profile considered) 

2 2 

Then it goes through the throat to another universe. The self-force will still attract the particle to the throat 
and the velocity of the particle will decrease. The particle will finally stop and the process will repeat. The 
process is similar to the oscillations of a particle in the potential well. Due to the acceleration the particle will 
lose its energy for radiation and in the end it will be at rest at the throat. Therefore particles will concentrate 
at the throat, at the minimum of the self-energy. The intensity of this process is managed by the characteristic 
value of the self-force, e 2 /a 2 , which is smaller for macroscopic wormholes. Therefore the smaller wormhole's 
throat, the faster it will gather neighbouring charged particles. At the same time particles at the throat with 
the same sign of charge will repel each other by Coulomb force. The magnitude of the force is defined by the 
same parameter e 2 /a 2 . Thus an equilibrium configuration of the particles near the throat is likely to exist. We 
can not exploit this picture for uncharged particles. Usually, the self-force for massive uncharged particle has 
opposite sign comparing with electromagnetic charge. The magnitude of the force is much smaller and defined 
by parameter Gm 2 instead of e 2 in the electromagnetic case. With reference to our case it means that the 
uncharged massive particle might be repelled by the wormhole. 
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